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&*' Abstract 

s 

The phenomenon of resonant activation of a Brownian particle over a fluc- 
tuating barrier is revisited. We discuss the important distinctions between 

q \ barriers that can fluctuate among "up" and "down" configurations, and bar- 

riers that are always "up" but that can fluctuate among different heights. A 
resonance as a function of the barrier fluctuation rate is found in both cases, 

O ■ but the nature and physical description of these resonances is quite distinct. 

The nature of the resonances, the physical basis for the resonant behavior, 
and the importance of boundary conditions are discussed in some detail. We 
obtain analytic expressions for the escape time over the barrier that explicitly 
^f) ' capture the minima as a function of the barrier fluctuation rate, and show 

that our analytic results are in excellent agreement with numerical results. 
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I. INTRODUCTION 

Noise-induced non-equilibrium phenomena in nonlinear systems have recently attracted 
a great deal of attention in a variety of contexts |1[. In general, these phenomena involve 
a response of the system that is not only produced or enhanced by the presence of the 
noise, but that is optimized for certain values of the parameters of the noise. One example 
is the phenomenon of stochastic resonance ||, wherein the response of a nonlinear system 
to a signal is enhanced by the presence of noise and maximized for certain values of the 
noise parameters. Another involves nonequilibrium ratchets, wherein intrinsically unbiased 
Brownian motion in stochastic asymmetric potentials leads to a systematic drift motion 
whose magnitude and even direction can be tuned by the parameters of the noise [§,|J . A 
third is the recent discovery of a re-entrant noise-induced phase transition in a nonlinear 
coupled array, that is, a transition that is only observed for certain finite ranges of noise 
parameters |J. A fourth such phenomenon, the one of interest to us in this paper, has been 
called "resonant activation" and was first identified by Doering and Gadoua and further 
studied by a number of other authors. Here the mean escape time of a particle driven by 
(usually white) noise over a barrier of randomly fluctuating height exhibits a minimum as a 
function of the parameters of the barrier fluctuations. 

Our interest in this problem first arose because it seemed to us that for sufficiently 
simple potentials it should be possible to find analytic dependences of the escape rate on the 
system parameters (or at least good approximations to them) and, more specifically, that 
it should be possible to find analytic expressions for the parameter combinations that lead 
to the minimum in the escape rate. Some analytic results are available |6|-|ll|, including 
those in the original work of Doering and Gadoua that apply to a very specific circumstance 
discussed in more detail below. In general, however, most available results are numerical 



I T[ . Analytic results are scarce, and usually apply only to one parameter regime or another 
and are thus unable to explicitly capture the occurrence of a minimum in the escape rate. 
A minimum in the escape rate usually arises from these approximations only by inference, 
and the approximations provide no way to specifically locate the minimum except as an 
intersection point of two unrelated approximations. They also do not provide a way to 
determine the dependence of the minimum on the system parameters. 

In this paper we accomplish our goal, that is, we obtain a number of analytic results for 
moments of the first passage time over a fluctuating barrier for the particular model system (a 
triangular potential barrier subject to dichotomous fluctuations) used in a number of studies 
of resonant activation. In particular, we obtain analytic approximations that explicitly 
capture not only the minimum in the escape rate but that allow us to study the variability 
of the escape rate in parameter space, that is, the depth and width of this minimum. 

In the process of obtaining these results we have also accomplished a number of important 
clarifications on the nature of models that have been presented under the common "resonant 
activation" rubric, and on the nature of resonant activation itself. Some of these models 
in fact differ from one another in essential respects. We discuss these clarifications and 
differences in some detail, and thus shed some light on the role played by the interplay of 
the white noise and the barrier fluctuations on the escape process. We anticipate some of 
our findings: 

• A distinction must be made between situations in which the fluctuating barrier can 



be "up" or "down" (i.e. can go from being a barrier to being flat or even a well), 
and situations in which there is always a barrier. Although resonant behavior can 
be observed in all cases, the physical picture underlying this behavior is different in 
different cases. 

• Boundary conditions play an extremely important role in the problem. 

• The qualitative physical description of the resonance in the fluctuating barrier problem 
is as follows. When the barrier fluctuates extremely slowly within a range < v m i n < 
v < v max , the mean first passage time to the top of the barrier is extremely long 
because it is dominated by those realizations for which the barrier starts in the high 
position. The mean first passage time is then proportional to e Vmax ' D where D is the 
intensity of the white noise. Indeed, if the barrier fluctuation rate is smaller than the 
inverse of the mean first passage time to the highest barrier, the barrier is essentially 
quasi-static throughout the process. At the other extreme, if the barrier fluctuates 
very rapidly, the mean first passage time is determined by the average barrier vq, i.e. 
it is proportional to e v °^ D . Between these extremes, and over a broad range of barrier 
fluctuation rates, passage over the barrier occurs primarily when the barrier is low, and 
the mean first passage time is then proportional to e Vrnin ' D . This dependence is quite 
robust, and the prefactor determines the actual minimum within this broad range. 

• This behavior does not require that the barrier fluctuate; an oscillatory variation of 
the barrier height yields essentially the same results. 

In Sec. |IT] we provide a detailed statement of the resonance activation problem. Section [D I 
discusses the analytic solution of the "up"- "down" case; we show that the resonance flipping 
rate and the resonance activation in this case are independent of the white noise intensity. 
In Sec. |V| the significance of the white noise and of the boundary conditions in this "up"- 
"down" problem are discussed in detail. Section deals with the case of barrier fluctuations 
when the barriers are always high. We obtain a single analytic expression for the mean first 
passage time that exhibits a minimum as a function of the barrier fluctuation rate and 
that in fact quantitatively captures the correct behavior over most of parameter space, as 
determined by comparison with numerical results. With this result we are able to determine 
the resonance frequency analytically, and also the range of barrier fluctuation rates over 



which the mean first passage time is essentially flat. In Sec. [VI] we discuss the case of a 
barrier that oscillates (rather than fluctuates). This case also exhibits resonant activation, 
although some of the quantitative details of the problem are slightly modified. Finally, we 
conclude with a summary and some final points in Sec. |V11| . 

II. STATEMENT OF THE PROBLEM 

Consider a process that evolves in a bistable potential and is driven by weak Gaussian 
white noise, so that the process is occasionally able to cross from one minimum of the 
bistable potential to the other. If the parameters of the system are fixed in time, the rate 
at which the process crosses from one well to the other under a variety of conditions is well 
known (e.g. the Kramers rate). Suppose now that the height of the barrier separating the 



two minima of the bistable potential fluctuates in time. We wish to explore the effect of the 
barrier fluctuations on the rate of passage of the process from one well to the other. More 
specifically, it is known that there is an optimal barrier fluctuation rate that minimizes the 
passage time from one well to the other for given parameter values |6|- |Tl|1 , This minimum 
identifies the phenomenon of resonant activation. We are interested in the analytic properties 
of the resonant activation phenomenon. Note that the barrier fluctuations here are such that 
the energy difference between the potential minima remains constant - - only the barrier 
height fluctuates. This is to be contrasted with the phenomenon of stochastic resonance, 
where the energy difference is modulated by a small periodic signal. 

We adhere to the overdamped regime, and hence the process y(r) evolves according to 
the Langevin equation 

y{T) = -V'{y)-g , (y)T}(T) + Z{r). (1) 

Here £(r) is zero-centered Gaussian white noise with correlation function 

<£(t)£(t')> = 2ZW(t-t'). (2) 

One can think of the white noise as arising from a heat bath, in which case the diffusion 
coefficient D is proportional to the bath temperature T. Time is measured in units of the 
friction coefficient, which has been set to unity in Eq. (TI|). 

The potential V(y) is a bistable potential, typically with isoenergetic minima. Doering 
and Gadoua |J introduced the triangular potential shown in Fig. |T|. The potential barrier 
is defined by 

v(v)-l Voy/L °<y< L (3] 

and the potential rises to infinity at y = and at y — 2L. In the absence of the contribution 
9'{y)v{ T ) i n Eq. (HD this represents a standard problem where the rate k at which the process 
crosses the barrier at x = L is related to the mean first passage time T\ from the bottom of 
one of the wells, say the one at y = 0, to the top of the barrier: k = \J2T\. To calculate 
the mean first passage time one assumes a reflecting boundary condition at y = and an 
absorbing boundary condition at y = L. 

In the resonant activation problem we have, in addition, the contribution g'{y)t](T). Here 
rj(r) is a non-equilibrium noise that coupled to g'(y) causes the potential barrier to fluctuate. 
It is a nonequilibrium noise because there is no dissipative contribution in the equation of 
motion associated with this fluctuating term, and hence the system is open. The noise 77 (r) 
is usually taken to be exponentially correlated, the most ubiquitous choices being Ornstein- 
Uhlenbeck noise PHTT| and Markovian dichotomous noise @JTT| . Here we deal only with the 
latter: tj(t) takes on the values ±1, and the change from one to the other is distributed in 
time according to the exponential density function 

0(r) = 7 e"^ (4) 

so that the flipping rate of the dichotomous noise is 7. The fluctuating barrier is accom- 
plished by picking for g(y) the function 



9 ^ = {-ay/l L +2a Liy<2L ® 

and zero otherwise. The addition of the random potential term g(y)r)(r) causes the potential 
barrier to switch between the two values v + = vq + ct and V- = vq — a. 

We wish to calculate the rate at which the process crosses the point y = L, which in turn 
is related, as before, to the mean first passage time Ti from y = to y = L when a reflecting 
boundary is located at y = and an absorbing boundary at y = L. (The distortions in the 
potential profile that may be caused by multiplicative noise, and the implications on the 
appropriate definition of an escape time, are well known and have been widely discussed in 
the literature; see. e.g. ITTTJ. The potentials used here do not exhibit such distortions.) In 



particular, we wish to establish analytically the dependence of T\ on the flipping rate 7 and 
to identify the flipping rate for which T\ is a minimum. 

Doering and Gadoua || calculated the mean first passage time for this model in the 
absence of the potential V(y), that is, when the "barrier" flips between being a true barrier 
(of height -\-a) to being a well (—a), and they obtained a resonance phenomenon, that is, 
the mean first passage time from to L exhibits a minimum at a particular value of the 
flipping rate 7. Doering and Gadoua also presented simulation results for the case v = a, 
that is, when the "barrier" flips between being a true barrier (height 2a) and there being 
no barrier. Bier and Astumian J7| considered the true barrier case, that is, the case where 
there is always a barrier (in fact, they took vq ^> a) and obtained numerical results that 
show a resonance. Their analytic barrier crossing rate results are obtained separately for low 
flipping rates (small 7) and for high flipping rates (large 7) . Neither result in itself exhibits 
a minimum, although one can infer the existence of a minimum (but not its dependence on 
the system parameters; see also |11|| ) from their combination. 

With this general statement of the problem we can be more precise about the results 
that we present in this paper. First, we consider the ±a barrier-well case of Doering and 
Gadoua, reproduce their analytic results for the mean first passage time, and also obtain 
analytic results for the resonant mean first passage time, the resonant flipping rate, and for 
the second moment of the first passage time distribution. We argue that the ±a barrier case 
represents a situation that is completely different from the "true barrier" case considered 
by Bier and Astumian. Both exhibit resonance behavior, but via different mechanisms. We 
explore these differences and interpret the Doering-Gadoua case on the basis of an even 
simpler model. Furthermore, we obtain analytic results for the high barrier case considered 
by Bier and Astumian that yield an explicit minimum in the mean first passage time as a 
function of the flipping rate. We present analytic results for the mean first passage time at 
resonance, for the resonant flipping rate, and we analyze the behavior of the system away 
from this point to assess how sharp this resonance might be. 

We shall present our analysis and results in terms of the dimensionless variables 

t = rD/L 2 , x = y/L, T^f.D/L 2 , (6) 

and the dimensionless parameters 

a = a/D, V = v /D, X = -fL 2 /D, (7) 

The differential equation for the mean first passage time in all cases considered in this 
paper is given by 



dx A dx 3 

with the boundary conditions 



dATl - 2V n ^ + (V* - a 2 - 2A)^ + 2XV dTl 



dx 2 



dx 



2A 



dT x 



dx 



x=0 



(8) 



(9) 



at the reflecting boundary and 



d 2 Ti 



dx 2 



x=0 



T 1 (x = l)=0, 



(10) 



(11) 



- — 2V - — I (V 2 2\ j^i 

<ix 3 <ix 2 i! dx 



V'o 



(12) 



x=l 



at the absorbing boundary. 

A brief description of how this equation and boundary conditions arise is given in the 
Appendix. 

III. ANALYTICAL SOLUTION OF THE DOERING GADOUA MODEL 

Consider the mean first passage time to x — 1 (y — L) when the mean barrier height 
is Vo = 0. The solution as a function of the initial position x can in this case be given 
analytically: 



T 1 {x) = {x-1) 



2Xa 2 \i — sinh(/i) 



-■^(x + 1) 



2a 2 sinh[fi(x - l)/2] 
/i 4 a 2 + 2A cosh(/x) 



/x 3 a 2 + 2 A cosh (/i) /i 2 
x {a 2 sinh[/i(x + l)/2] + 2Asinh[/i(x - l)/2] + 2/xAcosh[/i(x + l)/2]} 
where we have introduced the symbol 



//= Va 2 + 2A. 



(13) 



(14) 



This result has been previously reported for the particular initial value x = [H]. We 
know from Ref. that T\(x — 0) = T\ exhibits a resonance with respect to A; we wish to 
establish the resonance flipping rate \ res and the behavior of the mean first passage time at 
this resonance point. 

The expression for the mean first passage time simplifies considerably when the dimen- 
sionless quantity fi is large (fj, ^> 1), which is the physically interesting weak-white-noise 
regime. Indeed, the only way that \x could be small is if a and A are small, that is, if (in 
dimensioned units) the white noise intensity D is greater than the barrier height a and 
greater than 27L 2 . The first condition renders the problem uninteresting - if barriers are on 
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average lower than the noise then one has an essentially free diffusion problem. The second 
condition requires a small system with a low flipping rate, again a very specific situation that 
is not particularly interesting in this context. The customarily interesting physical situation 
occurs when the white noise is weak compared to the barrier height, that is, when a> 1 
and this in turn leads to /i ~^> 1. We use these two statements of the "interesting regime" 
interchangeably. 

When /i is sufficiently large the following approximations are valid: 

cosh/x — 1 ~ cosh/i ~ sinh/i — ji ~ sinh/z ~ ~e M . (f5) 

If in addition 

A > aV a (16) 

then the result simplifies even further, and one finally obtains the following much simpler 
approximate expression: 

aV -2\e~n , a 2 (2-x) A(l-x 2 ) 2 

Tl(x) ~ 2A(a 2 + 2A) 2 + p + 2A)3/ 2 + ^T^T + ° {a & ] ' (1?) 

It can be shown that this expression as a function of A has a minimum at a finite value X res 
that to leading order is of the form X res ~ a. This minimum is identified among the roots 
of dTi(x)/dx = as the one that coincides with the minimum of the complete expression 
of Ti(x), Eq. ([03), as a — > oo. Explicitly, following || we set x = to further simplify the 
analysis: 

(a 2 -2A)a 2 2a 2 A „, 2 „, 

The extrema of T\ as a function of A obey the equation 

%{a 2 + 2A) - ^(a 2 - 2A) - \{a 2 + 2A) - % (a 2 + 2A) 1/2 = 0. (19) 

This equation can be solved perturbatively by taking A ~ a (Ao + \ia~ l + A2a~ 2 ), an expan- 
sion consistent with the fact that a ^> 1. Once this expansion is substituted in Eq. (|19"D, the 
following result is obtained for the resonant flipping rate as a function of a: 

a„,~-L + (i + -L ) + ! ( 3 + -L,I, a>>1 . (20) 

Note that the condition A ^> a?e~ a is thus satisfied. The value of T\ at the resonant flipping 
rate is 



T 1 (A = A res )^T res ~(2 + v / 2)^-(4 + 3v / 2)^+^ + -^)^ 



(21) 



provided a is large. For small values of A (A <^i a 2 e a ) one readily finds from Eq. ( ]T3|) with 
x = that 



when a ^> 1. 

In Fig. ^ the exact expression of T\ is compared with the approximation (|22"D when 
A — ► and with the approximation ( Jl8| ) when A — > oo. 

The second order moment of the first passage time distribution, T 2 , also exhibits a 
resonance. However, the resonant frequency of T 2 does not coincide with that of T\. This 
means that there does not exist a unique resonant frequency or universal scaling associated 
with the first passage time distribution. T 2 can be calculated in a similar way (albeit even 
more expansively) than T\. However, the full expression for T 2 (x) is too long to be included 
here. Instead, we only reproduce the expression for T 2 at x = when a ~^> 1, 

m a 4 (a 4 + a 2 A - 4A 2 ) 4a 2 (a 2 - A 2 ) 2a 2 (4a 2 - A) 20Aa 2 5A 2 „, , u , , N 

T 2 ^ r^ - + — V^ — - + — — R + -5-r- + w~a + ° a e )• 23 

A 2 // 8 A/i 7 /i 6 3/x 5 3/x 4 v / v / 

In Fig. [5] this approximate result is compared with the exact expression for T 2 when a = 8. 
The behavior of the resonant frequency A reSi2 of T 2 can be obtained from the expression for 
T 2 when a — > 00. The calculation yields 

A res ,2~&oa + &i + 0(a~ 1 ), (24) 

where 60 is the positive real solution of the equation 5a 4 , + 106q — 660 H — 3 = 0. Numerical 
solution of this equation leads to the value 60 — 0.825724 . . .. With this value, the coefficient 
61 can also be evaluated numerically and one obtains 61 = 3.56057. Note that A reSj2 > A r 
Finally, T 2 at the resonant frequency behaves as 1/a 2 when a —>■ 00: 



^res- 



T 2 (A r . e ,, 2 )~^ + ^ + 0(a- 4 ), (25) 

a z a d 

with c = 20.9521 . . . and a = -104.244 . . .. 

Equations ([Rf), (pOQ, (pl|), and ( p2|) are the main results of this section - they give the 
mean first passage time as a function of the barrier height (or well depth) a, provided a <^C e a , 
for all values of A. The resonant behavior of the mean first passage time as a function of 
the flipping rate of the barrier is clear and dramatic. It is useful to exhibit explicitly the 
limiting results in the original units so that the dependence on system parameters is clear. 
When the flipping rate is large (A ^> a 2 or 7 ^> a 2 /DL 2 ) the mean first passage time from 
y = to y = L grows as 



T\ — > T w hite — 777^ — 1 (26) 



as seen from Eq. fllTf ). Here 



2A:// 



a 2 



D eff = D + —. (27) 

This is the well-known result for the mean first passage time from to L for a freely diffusing 
particle with diffusion constant -D e //- In this limit the flipping barrier behaves simply as an 
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additional source of white noise of intensity ifjji- At the other extreme, when the flipping 
rate is very low, A — > 0, the barrier never flips as the process moves from to L. If the barrier 
is initially "down" (77(0) = —1), then it remains down and the process is simply diffusively 
driven toward the absorbing barrier by a constant force; as a/D — > oo, the motion of the 
system becomes increasingly ballistic. If the barrier is initially "up" (rj(0) = 1), on the other 
hand, it remains up and the process moves between the reflective barrier at y = and the 
absorbing barrier at y = L against a constant opposing force. The mean first passage time 
for such a process grows exponentially with the barrier height as D 2 e a l D ja 2 ||12| . In our 
calculations either initial configuration is equally likely. The average of these two possibilities 
is the leading term in Eq. (^) (the ballistic contribution is negligible): 

Between these two limits lies a regime in which the mean first passage time is a minimum. 
The minimum value occurs at the resonant flipping rate whose leading term for large values 
of a/D is 

(29) 



"" x/2£2 
The leading contribution to the mean first passage time at this flipping rate is 

L 2 

T res ~(2 + v / 2)— (30) 

a 

and thus decreases with increasing a. This result has the L dependence of a diffusive 
process, but the effective diffusion coefficient here is a and not D. Note that the product 

ires J- res ^ v -*- / • 

IV. A SIMPLER MODEL SHOWING RESONANT ACTIVATION 

A surprising observation about the results of the Doering-Gadoua model is that the 
resonant frequency 7 res and the mean first passage time at this resonant frequency, T res , 
do not depend on the white noise intensity D to leading order in a. This means that as 
a/D — ► oo, the resonant properties become increasingly independent of the white noise 
intensity. Indeed, the resonance therefore appears unaffected by and unrelated to the white 
noise; in particular, the minimum in the mean-first-passage time in the Doering-Gadoua 
model appears not to arise from the coupling usually invoked between the white noise and 
the random dynamics of the potential. 

We have used somewhat equivocal language in this description because the situation is in 
fact somewhat subtle. The resonance in the Doering-Gadoua model arises from two features: 
1) the random dynamics of the potential (that is, the random switching between barrier up 
and barrier down) and, in particular, the initial average over these dynamics, and 2) the 
nature of the reflecting boundary at x — 0. It is this latter feature, subtly influenced by 
the white noise, that is especially noteworthy: the resonance characteristics of the Doering- 
Gadoua result when D — > are not reproduced by simply setting D = to begin with in 
the model equations. 
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To pursue this issue in more detail, let us consider the same model equations as did 
Doering and Gadoua but now in the absence of white noise from the outset. In place of 
Eq. ([ID, the system thus evolves according to the simpler Langevin equation 

y(r) = -g'(y)r]{r). (31) 

The solution of this mean first passage time problem is most easily found by splitting T(y) 
into two components |13|| : T + (y), the mean first passage time to y = L when 7/(0) = +1, 
and T~(y), the mean first passage time to y = L when ry(0) = —1. The convenience of this 
representation lies in the ease of expression of the boundary conditions in terms of T^. 

The boundary condition (H) is completely equivalent to the so called "immediate rein- 
jection" condition [[T|.|K|] 

f + (y = 0)=f-(y = 0). (32) 

In this case, whenever the system reaches the boundary at y = the velocity immediately 
changes its sign, that is, the driving noise t](t) changes its value from —1 to +1. Note 
that here the boundary condition directly affects the dynamics of the dichotomous barrier 
fluctuations since arrival at the boundary causes the noise to change its value. It is easy 
to ascertain that in terms of T(y) = (T + + T~)/2, the "immediate reinjection" reflecting 
boundary condition indeed translates to the Doering-Gadoua condition (P), i-e., 



df(y) 



dy 



= 0. (33) 



With this boundary condition (together with the absorbing condition at y — L) the mean 
first passage time from y = to y = L if the initial values 7/(0) = ±1 are equally probable 

is 



o 



Ti,ir = ^r + — • (34) 

a 2 a 

Note that T\^ r is a monotonically increasing function of 7 and thus exhibits no resonance. 
Clearly, this solution is not the one approached by the Doering-Gadoua model when D is 
set to zero in the solution of the latter. 

There is another way to think of a reflecting boundary, namely, to assume that the 
boundary only limits the region of movement of the system without interfering with the 
dynamics of the dichotomous barrier fluctuations. The dichotomous noise evolves according 
to its own dynamics, and changes its value at random times that are independent of where 
the process y{r) happens to be. Thus, if the system reaches the boundary y = when the 
noise happens to be 77 = —1, the noise may retain this value according to its own statistical 
properties. The process simply waits at the boundary, until the noise switches to 77 = 1 
in the natural course of events. We call this condition a "natural" reflecting boundary 
condition. This behavior is implemented via the following boundary condition for the mean 
first passage time components: 

f-(y = 0)=f + (y = 0) + - (35) 

7 
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or, in terms of T(y), 



dT(y) 



dy 

The solution for the mean first passage time is now 



„=-§■ (36) 

j/=0 



7L 4 2L 2 1 
a z a 27 



^l.n = V + + — • (37) 



It is easily seen that T\ >n has a minimum at 7 res = a/ \/2L 2 (see Eq. (|29|) ). In Fig. |] we 
have plotted a realization of the process y{j) for the two reflecting boundary conditions, the 
"immediate reinjection" and "natural". From this figure, it is clear that the two boundary 
conditions lead to different results for the mean first passage time. 

The interesting point to note is that the D — > limit of the mean first passage time in 
the Doering-Gadoua model is Ti >n , that of the "natural boundary," and not Ti iir , although 
the reflecting boundary condition used for the solution of the Doering-Gadoua model is 
(§). In the Doering-Gadoua model, no matter how weak the white noise, its effects become 
dominant near the reflecting boundary. The white noise allows reversal of the trajectory even 
infinitesimally close to the reflecting boundary without triggering a flip of the barrier, thus 
removing the dynamical interference between the boundary and the dichotomous flipping 
process. There is therefore a profound difference between the situation in which white noise 
is present, albeit as weak as one wishes, and the situation in which there is no white noise 
in the first place — the limit D — > is discontinuous. In the Doering-Gadoua model there 
is a resonance in the mean first passage time to absorption because, through the action of 
the white noise (no matter how weak) the process can delay (increasingly as the white noise 
intensity decreases) arrival at the absorbing state if the barrier is up and does not flip of its 
own dynamics. However, the characteristics of the resonant flipping rate and the mean first 
passage time at resonance do not explicitly depend on the intensity of the white noise. 

In Fig. |5] we exhibit the two mean first passage times, T\^ r and Ti n . 

V. THE BIER-ASTUMIAN MODEL 

The main conclusion that follows from the discussion of the preceding sections is that 
the resonant effect in the "toy" model of Doering and Gadoua is not of the same nature 
as the resonant activation in systems where the activation process is exclusively due to the 
presence of white noise (i.e., nonzero temperature). In order to study the resonant process 
in this latter situation, we return to the full model introduced by Doering and Gadoua but 
now with Vq > a, so that there is always a barrier. This problem was first analytically 
studied by Bier and Astumian [7J. The approximation developed by these authors coincides 
with the so called kinetic approximation introduced in Ref . [|TJ| . The main limitation of this 
method for the present purposes is that it leads to a mean first passage time that does not 
exhibit a minimum. 

We have developed an approximation for the mean first passage time to the absorbing 
boundary for high average barriers (Vo 3> 1) that does lead to a minimum and hence 
can be used to describe the resonance phenomenon analytically. We return to the Langevin 
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equation (|l|) with (0)-(|5]). Now, however, we take vo±a ~^> D or, in dimensionless quantities, 
Vo ± a ^> 1. The barrier thus flips between two large values. 
The general solution to Eq.(|8|) is 



T(x) 



Vo 



[x - 1) + A^e^ -e gi ) + A 2 (e q2X - e q2 ) + A 3 (e q3X - e 



q-.i \ 



where the coefficients Oj are the three roots of the polynomial equation 



q 3 - 2V q 2 + (V 2 



2X)q + 2AV = 



and the constants A{ have to be found from the boundary conditions 



)-m 



(38) 



(39) 
It can be 



demonstrated that for Vq > a the roots of Eq. (39) are all real, two of them positive and the 
other one negative. The full expressions for the constants Ai are complicated and too long 
to be included here. However, it is possible to derive shorter useful expressions for them as 
a series in A. In this case, the roots q^ can be written as 






2V 



ViV. 



A + 



-4 



4V (a 2 + V 2 ) 2 



2a Vi 



A VI 

v+ + vl + 2^V! 



A 2 + 0(A 3 ), 
A 2 + 0(A 3 ), 



where 



V± = V ± a. 



(40) 



(41) 



When these expressions are introduced into Eq. 
mean first passage time to order A 2 : 



T(x = 0) 



]), the following result is obtained for the 

N ± e v - + N 2 e v + + N 3 e 2V ° 



D x + D 2 e v - + D 3 e v - 
where the coefficients in the numerator are 

'VIV, 



(42) 



Nx 

N 2 
N 3 



Vl 



Yi 



A 
A 



'vv 



V 2 
aV_ 

V 2 



3aV+ a 2 V+ a(V -5a)' 



4A {1 + A 



V aV + 
(2V - 1) 

v+v_ 



v_ 

3aV_ 

4a 2 
(V+V-) 



v v_ vi 

a 2 V_ a( V) + 5a) " 



V V+ 



+ 



vl 



(43) 



and those of the denominator are 

Dx = 2(V+V_) 2 + 2A(a 2 + 3V - 2V V+V 
Do 2\ - ' 



Vl + 



(aV*V_ 



( ,,,.,_,_,_ V+V 2 + 2a 3 - 6a 2 V 
av_ 



D, = 2A 



V 2 + 



A 



aV 2 



aV 2 V + + V_V 2 + 2a 3 + 6a 2 V 



(44) 
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This approximation is one order higher in A than the one derived in Ref. |7j, which is 



equivalent to the so-called kinetic approximation |Tj| . The advantage of our new approxima- 
tion is that it shows a minimum as function of the frequency A. When V$ ^> 1 the resonant 
frequency can be calculated explicitly: 



A, 



a{e 2a - 1) 2 V 3 



4e 2a (l + a - e 2a + ae 



2« > 



1/2 



exp 



and the associated minimal mean first passage time reads 

2 



rp 

1 res ~ Tr2 i „-2aT/2 



V 2 + e~ 2a Vl 






(45) 



(46) 



The analytic expressions (^)-(|4q) are the principal results of this paper. 
As we did in the Doering-Gadoua case, it is useful to exhibit explicitly various limiting 
results in the original units so that the dependence on system parameters is clarified. 
When the flipping rate is large (A ^> a 2 or 7 ^> a 2 / DL 2 ) the result (|^) reduces to 



2 + a 2 /\ ( 2V 

1 1 —* 1 white — ^779 ex P 



2V 2 



2 + T- 



(47) 



or, in the original units, 



TWT, 



L 2 D, 



white 



v 2 o 



ff vo/D e ff 



(48) 



where -D e // is the effective diffusion coefficient defined in Eq. (p7|). This is the appropriate 
and familiar result for activation over a barrier of height vo with diffusion coefficient D e ff. 
At the other extreme, as A becomes small, the kinetic approximation [0,0] is valid and the 
mean first passage time (f42|) reduces to 



T^Tt 



km 



2X + (k + + k_)/2 
k + k_ + \(k + + k_) 



(49) 



where 



ku 



V 2 e~ v ±. 



(50) 



If A becomes so small that the time scale of barrier fluctuations is much slower than the 
escape time, then this further simplifies to 



TWT, 



hin 




(51) 



which is just the arithmetic mean associated with the two possible initial barrier configura- 
tions (see discussion surrounding Eq. (p8|)). In the original units 



Tl 



■ kin 



L2D e v + ,n + L 2 D 



,v-/D 



2v\ 



2vt 



(52) 
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Between these two limits lies the resonance regime where the mean first passage time is 
shorter than either the "white" or the "static" results. In the original units the mean first 
passage time at resonance, Eq. ([46]), reads 



f ~ 2L * D e"-"> 

res (vl + e-^vl) 

~ 2 -^-e v -/ D (53) 

vi 

where the second line, valid if a > 1, serves to stress the point that the resonant mean first 
passage time is essentially the usual passage time over the lower of the two barriers. It is 
not particularly instructive to exhibit the full expression fl45|) for the resonance frequency in 
the original units, but if a > 1 we can display the shorter expression 

3/2 

~ ~ —^ e ~ v - /2D (54) 

7res 2L 2 D 1 / 2 ' 

It should be noted that both the resonant mean first passage and the resonant frequency 
depend on the intensity of the white noise, as does their product. This dependence appears 
in the exponents as well as prefactors. 

A general feature of our solution and, more generally, of the resonant activation phe- 
nomenon is that with increasing barrier height the resonance phenomenon becomes less 
and less sharp: a long flat region develops around the resonant frequency, a fact that has 
been explicitly noted in earlier work ||17|| . Analysis of Eq.(fE^) makes it possible to estimate 
analytic bounds of this flat region, which spans the range 

V 2 e -2a + y2 yy_ 

—a e <<A «W ^ 

Thus, rather than stressing the resonance aspect of the problem, it might be more accurate 
to describe the time scale of the activation process as relatively insensitive to the parameters 
of the system except in the limits of very low and very high barrier fluctuation rates. As 
noted above, if the barrier fluctuations are sufficiently slow, then an initially high barrier 
remains that way essentially forever and the system on average crosses it before the barrier 
flips. Passage over the higher barrier then dominates the mean first passage time. At the 
other extreme, when the barrier fluctuations are very rapid, crossing occurs essentially over 
the average barrier. However, over most parameter ranges the mean first passage time is 
essentially determined by passage over the lower barrier - the system can avoid passage over 
the higher barrier by "waiting" for it to flip. Provided the waiting time is shorter than the 
time it would take the system to cross the high barrier, flipping will occur first and the 
system will cross when the barrier is lower (unless flipping is too rapid). This process is 
most efficient (but not dramatically more efficient - hence the flat behavior) at the resonance 
frequency. 

In Fig. ^ the mean first passage time and the different approximations explained above 
have been plotted for Vo = 11 and a — 1. Our approximation clearly captures the resonance 
behavior extremely accurately and for that matter the behavior of the mean first passage 
time over a broad range of barrier fluctuation rates. Fig. [7] again shows the mean first passage 
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time, but now with Vq = 15 and a = 1. This figure serves to reconfirm the agreement of our 
results with the numerical ones, and also illustrates the flattening of the region around the 
resonance as the height of the barrier increases. Fig. [8] shows the resonance frequency as a 
function of the barrier height. Again, our approximation clearly captures the exact results 
extremely well for sufficiently high barriers. 

VI. ACTIVATION DRIVEN BY A SQUARE WAVE FUNCTION 

It is interesting to explore whether the resonant activation phenomenon requires that 
the barrier fluctuate stochastically, or whether it also occurs when a noisy process crosses 
a barrier that changes periodically. Indeed, stochastic fluctuation of the barrier is not a 
requirement. 

To investigate the activation process when the barrier oscillates periodically between a 
higher value and a lower value, we replace the dichotomous noise r)(r) in Eq. ([I]) with a 
square wave w(t), a periodic function that alternately takes on the values +1 and —1. The 
changes from one to the other occur at a constant frequency 7. The period of the square 
wave function is thus 2/7. 

The Fokker-Planck equation describing the evolution of the probability for the system 
now includes a time-periodic potential. The problem can be tackled analytically using 
Floquet theory. We simply state qualitatively the results that one obtains with this exact 
approach, but then follow a simpler approach to arrive at some quantitative conclusions. 

We continue our discussion in terms of dimensionless variables and parameters. Exact 
solution of the problem does, as noted above, also lead to resonant activation when the 
barrier changes from higher to lower periodically, with a resonance behavior very similar 
to that of the stochastic case. In other words, the mean first passage time is large when 
the period of oscillation is very slow and also when it is very fast. As before, and for the 
same physical reasons, in the former case the mean first passage time is dominated by the 
high barrier V + and in the latter case it is determined by the average barrier Vo- Again as 
before, between these two limits there is a flat region (i.e., rather insensitive to the parameter 
values) where the mean first passage time is determined primarily by the lower barrier VL. 
The only difference between this problem and the stochastic one lies in the detailed way in 
which the mean first passage time changes from one behavior to the other. 

To find the mean first passage time at the slow-barrier-modulation end of the problem 
(where the difference between stochastic and periodic modulation is most pronounced), we 
recall that for a fixed barrier of height V the probability that the process has not yet crossed 
the barrier at time t (i.e., the survival probability at time £), is exponential [H|, e~~ kt , where 



the crossing rate k = V 2 e~ v [cf. Eq. (0)]. If the barrier is not fixed, but instead changes 
slowly from one value to the other at predetermined periodic time intervals A = A -1 with 
A < 1, we can track the trajectory of the system explicitly and write down an expression for 
the survival probability S(t) that the process has not yet crossed the barrier at time t. If 
initially the barrier is V+ and the time t = corresponds to the beginning of a period, then 

S + = Prob{T+ >t} = qlq^e- k+{t - 2nA \ 2nA < t < (2n + 1)A; 

S+ = Prob{T + >t} = g ™+ig™ e - MM2 " +1)A] , (2n + 1)A < t < (2n + 2) A (56) 
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where n = 0, 1,2, . . ., and q + and g_ are the probabilities that a crossing event does not 
happen when the barrier is, respectively, V + and VL. The assumption about the statistics 
of the crossing events yields 

q ± = e -*±/A_ (57) 

The mean first passage time T + can then be calculated directly as a moment of this probabil- 
ity. The survival probability S~ and associated mean first passage time T~ when the barrier 
is initially V_ is similarly obtained. To compare most directly with the stochastic results we 
assume that initially the barrier is equally likely to be V + or l/_. A short calculation then 
leads to the following result for the mean first passage time when A ^$ 1: 

T+ + T~ 1/1 1 \ 1/1 1 \ q+ - q- 

r =^^=o ir + ir +o r-rrrr- ( 58 ) 



This result corresponds to the same level of approximation as the kinetic result (|49|) . At 
very low frequencies, A — ^ 0, the mean escape time is correctly given by an average of the 
escape time l/k + when the barrier is V+ and 1/k- when the barrier is V-. This kinetic 
approximation also does not exhibit a minimum because it does not behave correctly when 
A — > oo; instead, it converges to the same value as the kinetic approximation ( fi9|) , that is, 
to (k + + fc_)/2. 

In Fig. ^| the approximation (|58|) to the mean escape time for the activation process 
driven by a periodic signal is compared with the escape time for the same system driven 
by dichotomous noise. The difference between the two is noticeable in the decrease of the 
mean first passage time with increasing flipping rate -- the dependence on flipping rate is 
considerably sharper in the periodic case than in the random case. A similar effect has been 



observed recently in systems that exhibit coherent stochastic resonance [|19|]. The minimum 



first passage time and resonance flipping rate are essentially identical in the two cases. 



VII. CONCLUSION 

We have revisited the problem of resonant activation, that is, of the mean escape time of 
a particle driven by white noise of intensity D over a barrier of randomly fluctuating height. 
The initial position of the particle is y = and the barrier is at y = L. A substantial recent 
literature [&-pi O , ^B[ deals with this problem, but the results to be found in the literature 
are almost exclusively numerical. The distribution of barrier fluctuations is typically taken 
to be either dichotomous (i.e., the barrier fluctuates between two values) or Gaussian. The 
correlation function of the barrier fluctuations is usually assumed to be exponential and thus 
characterized by a rate parameter 7. The quantity of interest is the mean escape time T\ 
of the particle over the barrier as a function of 7. It is observed that T\ vs. 7 exhibits a 
minimum, i.e., there is an optimal barrier fluctuation rate that minimizes the escape time 
of the particle. This minimum defines the resonant activation phenomenon. 

In this paper we have concentrated on dichotomous fluctuations and on triangular po- 
tential barriers, so our quantitative results are restricted to these cases. However, we believe 
that our results provide insights beyond these specific conditions. In particular, they pro- 
vide insights for barriers whose fluctuations are bounded between an upper value v + and 
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a lower value V-. The applicability of our conclusions to Gaussian fluctuations is therefore 
less certain, but below we will present some conjectures for this case as well. 

A variety of approaches to the problem of the escape over a fluctuating barrier of bounded 
variation provide excellent and consistent analytic approximations to the escape rate in the 
limiting cases of very slow barrier fluctuations and of very fast barrier fluctuations. In the 



slow fluctuation case, the so-called "kinetic approximation" [16| captures the behavior of 
the system very well. In the limit of very slow fluctuations (7 — ► 0) the barrier retains its 
initial height throughout the process. The mean first passage time for the ensemble is then 
just the mean first passage time averaged over the initial distribution of barrier heights. 
For example, in the dichotomous case if the height of the high barrier is v + then the mean 
escape time as 7 — > is determined by the mean escape time over this high barrier (the 
mean escape time over the lower barrier being negligible in comparison): 

f static ~ ^r^' D . (59) 

This is the result captured, for instance, in Eq. ( ffiD and in Eq. Q5"2| . Clearly, this result 
is determined in part by the assumption (generally made in the literature) that an initial 
average over an ensemble of barrier heights is appropriate. The entire discussion that follows, 
including the occurrence of a resonance, is dependent on such an initial average or at least 
on the assumption that a finite fraction of realizations begin with a barrier configuration 
that is higher than the lowest barrier. 

At the opposite extreme, when the barrier fluctuations are very rapid (7 — > 00), the 
main effect of the flipping barrier is to increase the effective intensity of the white noise. 
The escape then occurs over the average barrier, with a diffusion coefficient D e ff which 
exceeds D by an amount determined by the detailed distribution of barrier fluctuations. If 
the average barrier height is vq > then the mean escape time in this limit is 

f^ue-^^e^^ff. (60) 



vl 



If the average barrier height is zero, then 



T wh ite - 777; — • ( 61 ) 



These are the results captured in Eqs. (pq) and fl48|). In any case, the escape time is clearly 
smaller in the fast barrier fluctuation limit than in the slow barrier fluctuation limit. 

In the literature, each of the above approximations had been carried sufficiently far to 
deduce the behavior of the escape time as one moves away from the strict limits. Thus, within 
the kinetic approximation it can be shown that the escape time decreases with increasing 
7. At the opposite limit, it can be shown that the escape time decreases with decreasing 
7. These two results clearly point to a minimum for some finite value of 7, but neither 
approximation is sufficient to actually capture the minimum. Our goal here has been to 
develop a single approximation to capture this minimum, and in this we succeeded. 

However, we found in the process that a distinction needs to be made between two 
cases that lead to a different physical origin and parameter dependences for the resonant 
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flipping rate and the associated escape time. These two cases have in the literature been 
treated more or less as one because both involve dichotomous fluctuations, but they are in 
fact very different. One of these is the case in which the "barrier" fluctuates between an 
"up" or positive (barrier) configuration of height v + = a and a "down" or negative (valley) 
configuration of height v_ = —a. We have called this the Doering-Gadoua model ||. In 
the other case, the barrier fluctuates between a high value v + = vq + a and a lower (but still 
positive) value V- = vq — a. We have called this the Bier-Astumian model M. For each 
model we found a single expression for the mean first passage time that has a minimum, 
and we compared our results with exact ones obtained numerically. The agreement in both 
cases is excellent for almost the entire range of flipping rates, and in particular over a broad 
range surrounding the resonance. 

The distinctive aspect of the Doering-Gadoua model is the fact that part of the time the 
"barrier" is really a valley, so that the particle can essentially roll rather than climb towards 
L during these times. We found an explicit expression for the resonant flipping rate and the 
resonant mean first passage time in this case: 

a 

Ires ~ Wl , (62) 

L 2 

T res ~(2 + V2)—. (63) 

a 

The noteworthy fact about these results is that neither the resonant flipping rate nor the 
resonant escape time depend explicitly on the intensity D of the white noise. This fact 
seems not to have been noted before. We then went on to explore whether in fact this 
resonance is observed in a process defined by the Doering-Gadoua model with no white 
noise from the outset, and found that there is no resonance for such a model! We explained 
this apparent contradiction by noting a discontinuity in the D — > limit of the problem and 
by presenting a modified set of boundary conditions that does lead to a resonance (precisely 
the Doering-Gadoua resonance) in the absence of white noise. 

In order to obtain a result for the escape time in the Bier-Astumian model that cap- 
tures the resonant behavior we found that we had to retain terms in our solutions to one 
power higher in 7 than had been done previously (the lower orders yielded only the kinetic 
approximation) 0. With this, we identified the resonant frequency and escape times as 

3/2 

~ ~ —^ e~ v - /2D (64) 

< res 2L 2 D 1 / 2 

and (more complete results are found in Sec. [V]) 

f res ~ ^e«-/» (65) 

vi 

We noted that that in this case both the resonance frequency and the mean first passage 
time at the resonance frequency depend on the intensity of the white noise, as does their 
product. At resonance the escape over the barrier occurs primarily when the barrier is at its 
lowest. We also noted that the dependence on the flipping rate, especially for high barriers, 



is very flat: there is a broad range of flipping rates where passage over the barrier occurs 
primarily when the barrier is low. In this broad range of flipping rates the escape time over 
the high barrier is so long that the barrier is likely to flip to its lower height before the 
escape is completed. 

We also discussed the fact that the resonant activation phenomenon does not require a 
fluctuating barrier — it also occurs if the barrier oscillates periodically between the high and 
low values. The behavior of the escape time at low and high oscillation periods is the same 
as in the dichotomous fluctuation case, and at intermediate oscillation periods a resonance 
effect is also observed. 

Finally, we note that our analysis does not address the case of Gaussian barrier fluctua- 
tions, that is, of Ornstein-Uhlenbeck barrier fluctuations fll ]. The results for such barrier 



fluctuations with fixed variance [T(J should be similar to our results for dichotomous noise. In 
particular, the escape time for the model analogous to that of Doering and Gadoua (Vo = 0) 
will show a minimum even in the absence of white noise. 
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APPENDIX A: SURVIVAL PROBABILITY AND FIRST PASSAGE TIME 

MOMENTS 

The moments of the first passage time from x = to x = 1 can be obtained from the 
survival probability S(x, t) that the system evolving according to Eq. (]T|) (appropriately 
scaled to dimensionless variables) with a reflecting boundary at x = and an absorbing 
boundary at x — 1 has not left the interval (0, 1) at time t. This survival probability obeys 
the following partial differential equation (a detailed derivation and original references for 
the survival probability in an interval terminated by two absorbing boundaries are presented 
in Ref. f2§): 

£ 2 S + 2X£S = a 2 ^r, (Al) 

ox 1 

where £ is the differential operator 

_ d d d 2 

dt dx dx 2 

Note that Eq. (|A1|) is a second-order partial differential equation in the time variable and 
a fourth order partial differential equation in the state variable. Therefore, two initial 
conditions and four boundary conditions are needed to solve it. The initial conditions are 

5(x,0) = l, (A3) 

f| - (-» 

Uh t=o 
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For the absorbing trap at x = 1, the boundary conditions read 

S(l,t)=0, 



'c + v°-v 2 + a 2 )^- 

ox j ox 



-V 8(t), 



x=l 



and for the reflecting boundary at x = they are 

dS_ 
dx 



0. 



x=0 



d d 2 



S 



x=0 



0. 



(A5) 



(A6) 



(A7) 



(A8) 



\ dt dx 2 
The first-passage-time moments T n are related to the survival probability according to 

T n {x)=n e^Six^dt. (A9) 

Jo 

Clearly Tq(x) = 1 by normalization. T\(x) is the mean first passage time to 1 for a process 
that starts at X(0) = x; T^x) is the second moment of the distribution, so that the variance 
of the distribution of mean first passage times is a 2 = T<i — T 2 . 

Equations for the first passage time moments can be obtained by multiplying Eq. (|A1|) 
by t n ~ l and integrating over time by parts. The following recursive-differential equation is 
easily found: 



^D-^n — 2\L,]ji± n 



2 Ct l n 

dx 2 



9m 



where 

g n = n(2\ - 2C D )T n _ 1 - n(n - l)T n _ 2 , 
and Cd is the Fokker-Planck operator 



£ 



D 



V— — 
dx dx 2 



(A10) 



(All) 



(A12) 



with T_i =0. The boundary conditions can be obtained directly from those of the survival 
probability: 



1. TJ1) 







2 ( c ° - < + ^ " A f 



dT, 



-n 



n-l 



x=l 



dx 



+ V 6, 



n.l 



x=l 



3. 



4. 



dT n 
dx 

d 2 T n 

dx 2 



x=0 



x=0 



-nT n _!(0). 



(A13) 
(A14) 
(A15) 
(A16) 
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FIGURES 
FIG. 1. Schematic of the fluctuating potential barrier problem. 

FIG. 2. Mean first passage time as a function of barrier fluctuation rate for the Doering-Gadoua 
model. The barrier fluctuates between the up and down positions with slopes a and —a respectively. 
Solid curve with circles: exact mean first passage time obtained numerically. Solid curve with 



squares: our analytic result (|18|) . Dotted curve: the low-frequency approximation (22). The 
parameter a = 8. 

FIG. 3. Second moment of the first passage time distribution for the Doering-Gadoua model 
as a function of barrier fluctuation rate, with a = 8. Solid curve with circles: exact second moment 
obtained numerically. Solid curve with squares: our analytic result (p3j). 

FIG. 4. Typical trajectories near the reflecting boundary for the two types of boundaries dis- 
cussed in the text. Dashed lines: "immediate rejection" boundary. Solid lines: "natural" boundary. 

FIG. 5. Mean first passage time to L = 1 as a function of barrier fluctuation rate for the 
two types of reflecting boundaries discussed in the text. Dashed curve with squares: "immediate 
rejection" boundary. Solid curve with circles: "natural" boundary. 

FIG. 6. Mean first passage time as a function of barrier fluctuation rate for the Bier-Astumian 
model. The barrier fluctuates between the high value Vq + a and the lower value Vo — o> with 
Vq = II and a = 1. Solid curve with circles: exact mean first passage time obtained numerically. 
Dotted curve: kinetic approximation (ff9|). Dashed curve: white noise approximation (|47|) . Solid 



curve with squares: our result (42). 



FIG. 7. Mean first passage time as a function of barrier fluctuation rate for the Bier-Astumian 
model with Vq = 15 and a = 1. Solid curve with circles: exact mean first passage time obtained 
numerically. Dashed curve: white noise approximation (|47|). Solid curve with squares: our result 



FIG. 8. Resonance frequency as a function of average barrier height for the Bier-Astumian 



model with a = 1. Circles: exact numerical results. Solid curve: our result (45). 



FIG. 9. Mean first passage time as a function of barrier oscillation/fluctuation rate. Square 
symbols: kinetic approximation for square wave barrier oscillations. Circles: fluctuating barrier. 
Barrier parameter values: Vq = 11 and a = 1. 



22 








2L 



2.0 



1.5 



1.0 



o 



0.5 



0.0 



-0.5 




CM 



O 








LogioW 




Reflecting Boundary 



- 1 



o 








Logio(y) 





2.90 




2.85 




2.80 




2.75 


o 

D) 
O 

_l 


2.70 




2.65 




2.60 




2.55 




2.50 




Log 10 W 



4.30 



4.25 



4.20 



° 4.15 

O 



4.10 



4.05 



4.00 




eeduhhush 



-4 




LOQioW 



CO 



c< 



8 2 







5 



10 



15 



V, 



2.85 



2.80 



2.75 



° 2.70 

O 



2.65 



2.60 



2.55 




